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, We report a theoretical calculation explaining the quantum Nernst efEect observed experimentally 

^N) ' in a bismuth single crystal. Generalizing the edge-current picture in two dimensions, we show that 
the peaks of the Nernst coefficient survive in three dimensions due to a van Hove singularity. We 

QQ ^ also evaluate the phonon-drag effect on the Nernst coefficient numerically. Our result agrees with 

C " 3 1 the experimental result for a bismuth single crystal. 
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In the present study we extend the the- 
oretical prediction of the quantum Nernst 
effect in two dimensions -'^ to a three- 
dimensional system, motivated by the ex- 
periment in a bismuth single crystal. We 
find that the peaks of the Nernst coeffi- 
cient survive in three dimensions because of 
a van Hove singularity. A simple calculation 
explains Behnia et a/.'s experimental result 
qualitatively. Consideration of the clcctron- 
phonon interaction leads even to quantitative 
agreement with the experimental result. 

2. Edge-current picture of the 
quantum Nernst effect 

We consider the adiabatic Nernst effect. The 
magnetic field B is applied in the direction of 
the trigonal axis, which we set to the z axis. 
The system is electrically and thermally insu- 
lated on all surfaces except for the two sur- 
faces attached to heat baths, which induce 
a temparature bias AT in the x direction. 
Then the Nernst voltage Vn is generated in 
the y direction. The adiabatic Nernst coeffi- 
cient is then defined by 

N=-^, (1) 

where L and W are the lengths of the crystal 
in the x and y directions, respectively. 

Nakamura et al} assumed that, in a two- 
dimensional clean sample at low tempera- 
tures, an edge current circulates around the 
system ballistically when the chemical poten- 
tial is in between two Landau levels. The 
edge current that departs the colder side 
maintains the Fermi distribution /(T_,/u_) 
of the colder bath until it reaches the hot- 
ter side and then maintains the Fermi dis- 
tribution /(T+,/i+) of the hotter bath until 
it reaches the colder side. The difference of 
the chemical potential generates the Nernst 
voltage as Vjv = — /x_)/|e|. 

The Hamiltonian of the two-dimensional 
electron gas yields energy levels E{n,kx), 
where n = 0, 1,2, ••• represents a Landau 



level and the label — fcm < fcx < repre- 
sents a state in the Landau level with fcm de- 
termined by the confining potential in the y 
direction. After some algebra they obtained 
the electric current in the x direction as 

/ = 4 (^oA/z + AikeAT) , (2) 
nil 

where AT = — A/k = — /i_ and 
•Xi(n) ^^^^ 



o-Jxo(n) 4 cosh (x/2) 



(3) 



with Xi{n) = {Ei{n) - n)/{kBT), Eo{n) = 
E{n,kx = 0) and Ei{n) = E{n,kx = km)- 
Setting 7 = and e = |e| > 0, they obtained 
the Nernst coefHcient as^ 

ks L Ai 



N{B,T) 



\e\B,WAo' 



(4) 



which showed sharp peaks when the chemical 

potential is equal to a Landau level. This 
quantum behavior basically comes from the 
oscillation of the coefficient Ai . 

We now extend the above to three di- 
mensions. This extension introduces an ad- 
ditional argument kz in the energy level as 
Ei{n, kz) = Ei{n) + {hkz)'^ / {2mz) and hence 
in Xi{n,kz)- Equation (3) now has integra- 
tion over kz as 



CO 



M^^) = y:i^ 



n=Q 



2tt 



Xi(n,/c^) 



-f 

Jo 



Xa(n,k,) 4 cosh (x/2) 
D{Ez)A,{,i-Ez)dEz, (5) 



where we have changed the integration vari- 
able from kz to Ez = {hkz)'^ /{2mz). This 
variable transformation introduces the den- 
sity of states D{Ez) = [2'K{dEz/dkz)]~^, 
which has a van Hove singularity at = 0. 

This van Hove singularity is the essence 
of the survival of the quantum oscillation in 
three dimensions. The overlap of the peaks of 
with the van Hove singularity leads to 
oscillatory behavior of the coefficient 
and hence of the Nernst coefficient. 

Figure 1 presents the quantum Nernst ef- 
fect thus obtained in three dimensions. This 
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Fig. 1. The Nernst coefficient A'^ times the magnetic 
field, which is equal to the thermomagnetic power 



in three dimensions at T = 0.28 K against the 



inverse magnetic field 1/B [1/T], after Eq. (4). We 
used the following parameter values: the effective 
mass is m = 0.06289mo and rriz = 0.6667mo^ where 
mo is the bare mass of the electron; the size of the 
Hall bar is L = 4.0mm'''5 and W = 2.2mm-*'5; the 
confining potential at the edges is Vo = 4.22 eV'^; the 
chemical potential at equilibrium is ^ = 11.4meV^. 

agrees with the experimental result qualita- 
tively, though the magnitudes of the peaks is 
much smaller in our calculation than in the 
experimental result. 

3. Phonon-drag effect on the 
Nernst coefficients 

Since the above argument does not quanti- 
tatively agree with the experimental result, 
we now consider the phonon-drag effect. We 
basically follow Sugihara's theory^'^° of the 
thermoelectric power in bismuth. The differ- 
ence here is that we treat the Fermi and Bose 
distribution functions without any approxi- 
mations and evaluate the magnetic-field de- 
pendence of the thermomagnetic power nu- 
merically. 

The thermomagnetic power Sxy — NB 
is given by 

Qy 



^xy — rppPxj: 



(6) 



where Qy denotes the heat current in the y 
direction, F denotes the electric field in the x 
direction, p^x denotes the diagonal resistiv- 
ity, and we assumed pxy <C Pxx in bismuth. 
At low temperatures we may neglect all ex- 
citations except acoustic phonons having the 



energy fkOq and the wave vector q, which 
are generated through deformation coupling. 
The heat current in the y direction is then 
given by 



Qy 



dq 



fuVqVs — ( N, 



q 



(7) 



where Vs is the group velocity of the phonons. 
The term Nq — Nq^^ represents the displace- 
ment of the phonon distribution Nq from its 
equiribrium Bose distribution Nif "' . 

In order to estimate the displacement, 
we use the Boltzmann equation in the steady 
state 



fdNq 



dN„ 



dt 7 carrier V dt 



0. (8) 



The first term represents the change of the 
phonon distribution due to the interaction 
with carriers and the second term repre- 
sents that due to other interactions such as 
boundary scattering, phonon-phonon inter- 
action and impurity scattering. In the steady 
state these two terms are balanced. In the 
second term, we assume (9iVq/9i) relaxation — 
— {Nq — N!f'^)/T-[{q) in the relaxation-time 
approximation; the electron-phonon interac- 
tion replaces the phonon distribution, but 
the relaxation effects make the nonequilib- 
rium phonon distribution to the equilibrium 
one in time t^-. 

The theoretical result is shown in 
Fig. 2 (a) and the experimental result is 
shown in Fig. 2 (b) at three temperatures, 
respectively. In Fig. 2 (c) we compared 
the theoretical and experimental result at 
T=0.28[K]. The magnitude of the peaks of 
the theoretical result, after considering the 
phonon-drag effect, is consistent with the ex- 
perimental result. Note that there are no ad- 
justable fitting parameters. 

4. Summary 

We have shown that the quantum Nernst ef- 
fect survive in three dimensions due to a van 
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Fig. 2. The thermomagnetic power in three dimen- 
sions at T= 0.28 K, 0.46 K, 1.20 K against the in- 
verse magnetic field l/iJ[l/T]: (a) Our theoretical 
result with the phonon-drag eff'ect; (b) The exper- 
imental result by Behnia et a/.*'^; (c) Comparison 
of the results at T= 0.28 K. In the calculation we 
used the following parameter values: the deforma- 
tion potential D = 1.2eV^^; the electric resistivity 
pxx = 0.35f2cm*'^; the group velocity of the phonons 
2 X lO-^m/s; the bismuth density 9.75 X lO^kg/m^. 
The parameter values for the effective mass m and 
and the chemical potential at equilibrium fi are 
the same as in Fig. 1. 



Hove singularity. We have presented an ex- 
tension of the edge-current picture in two di- 
mensions to three dimensions, in order to ex- 
plain Behnia et a/.'s* experiment for a bis- 
muth single crystal. The naive extension re- 
produced the quantum oscillation of the ex- 
perimental result qualitatively. Furthermore, 
consideration of the phonon-drag effect led 
us to a theoretical result quantitatively con- 



sistent with the experimental result. The 
peaks are sharper than those in the exper- 
iment probably because we neglected impu- 
rity scattering of carriers. For simplicity, we 
also assumed the Fermi surface of bismuth is 
ellipsoid and considered only the hole contri- 
bution. Going beyond these approximations 
would be an interesting future problem. 
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